Magnetic dipole and electric quadrupole moments of the 229 Th nucleus 
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We calculate the A and B hyperfine constants for the low- lying states of 229 Th 3+ using a high- 
precision relativistic all-order approach. By combining these calculations with measurements of the 
6d and 5/ hyperfine constants [C. J. Campbell et al, Phys. Rev. Lett. 106, 223001 (2011)], we 
determine the magnetic dipole /j, = 0.360(7) fiB and the electric-quadrupole Q — 3.11(6) eb moments 
of the 229 Th nucleus. Our value for fj, is five times more accurate and is 22% smaller than the best 
previous value fi = 0.46(4) [ib [S. Gerstenkorn et al, J. Phys. (Paris) 35, 483 (1974)], while our 
value for Q is the same, but 2.5 times more accurate than the 2011 result. A systematic study of 
hyperfine structure in eight other monovalent atoms supports our claim of 2% level uncertainty for 
fi and Q in 229 Th+ 3 . 

PACS numbers: 31.30.Gs, 21.10.Ky, 27.90.+b, 31.15.ac 



I. INTRODUCTION 

The 229 Th nucleus has an unusually low first excita- 
tion energy of only several eV [l], Q , making the corre- 
sponding nuclear transition accessible with laser excita- 
tion. The transition is expected to be very narrow and 
well-isolated from effects of external fields and, therefore, 
presents a remarkable opportunity for the development 
of a nuclear clock Q . The physical implementation may 
employ either the stretched states within the 5/5/2 elec- 
tronic ground level of both nuclear ground and isomer 
manifolds of a single trapped ion [H, or the closed elec- 
tronic shell of Th 4+ in a UV-transparent crystal doped 
with a macroscopic number of 229 Th nuclei [5[. More- 
over, the transition frequency is expected to be as many 
as five orders of magnitude more sensitive to temporal 
variation of the fine structure constant and the dimen- 
sionless strong interaction parameter m q / Aqcd as com- 
pared to atomic transitions, making 229 Th one of the 
most attractive candidates for such studies @. Th 3+ 
was also suggested for study of parity violation 0] ■ 

For the future development of the ion-trap nuclear 
clock, it is essential to possess accurate values for 
both nuclear and electronic properties of 229 Th 3+ . 
The electronic level structure of the Fr-like Th 3+ 
ion, with a single valence electron above the closed 
[Rn] = [Xe]4/ 14 5c? 10 6s 2 6p 6 core, is conducive to high- 
accuracy atomic calculations d, . Laser-cooled Wigner 
crystals of 229 Th 3+ allow for high-precision spectroscopy 
[lol . fllT ] . By combining measured hyperfine B constants 
for the 5/5/2, 5/7/2, 6d 3 / 2 , and 6d 5 / 2 levels with atomic 
structure calculations of Ref. ||, the nuclear electric 
quadrupole moment Q = 3.11(16) eb was reported [Toj . 

From a broader perspective, precision measurements 
of nuclear moments by atomic spectroscopy are essential 
to improve our understanding of the atomic and nuclear 



structure. Such measurements may also lead to an emer- 
gence of novel physics involving, e.g., breaking of a fun- 
damental symmetry or a temporal variation of a funda- 
mental constant. The on-going proton radius puzzle fl2| 
underscores the significance of such studies. In this Let- 
ter, we report (1) a determination of the magnetic dipole 
and electric-quadrupole moments of the 229 Th nucleus at 
the 2% level of precision, and (2) a benchmark study of 
the accuracy of a theoretical method for the calculations 
of the electric quadrupole hyperfine matrix elements for 
future determination of electric-quadrupole moments in 
other systems. We also discuss the prospects for fur- 
ther accuracy improvements and possible extraction of 
the nuclear magnetic octupole moment. 

The atomic hyperfine constants A are proportional to 
the nuclear magnetic dipole moment /1. Therefore, we 
calculate (A/ fj,) th and determine /x as the ratio of the ex- 
perimental value for the hyperfine constants A expt and 
(A/fi) th . Similarly, the nuclear electric quadrupole mo- 
ment Q is determined as the ratio of the experimental 
value for the hyperfine constants B cxpt and the computed 
quantity (B/Q) th . 

The calculations are carried out using a relativistic 
linearized coupled-cluster (all-order) method including 
single, double, and partial triple (SDpT) excitations 
described in review [13j and references therein. This 
method was applied to the calculations of transition 
properties and polarizabilities of Th 3+ in Refs. [H,[Hj]- A 
brief description of the method is given in the supplemen- 
tary material [151 ] . This method has been applied to the 
calculation of h ype rfine constants in a number of mono- 
valent systems jl6l - tl^ . [l^ - [2lT ] . The resulting values for 
Ca + , Sr + , Hg + , Rb, Cs, Ba + , Fr, and Ra + are compared 
with experimental values in Table I of the supplementary 
material fl5| . The nuclear magnetic moments are taken 
from [22j unless noted otherwise. The results for the last 
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TABLE I: Comparison of the hyperfine constants A (in MHz) 
with experiment in 133 Cs Ej, H, 137 Ba+ [H^, 201 Hg+ 
H|, 210 Fr (/n = 4.38(5)/x s ) |Iap5(fl J^Ra 4- (fi = 0.878(4)/^?) 
and 213 Ra+ (/x = 0.613(2)/x s ) [31H33|]. Theory values for Cs 
and Fr are from 120]. Ba + and Ra + values are calculated in 
the present work. Lowest-order DF values are given in col- 
umn DF and relative correlation corrections are given in col- 
umn "CC" in %. "Rescaled from the measurement in 223 Ra. 
^Recommended value based on Ba + data. 
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four systems are given in Table H We find that Ba+ and 
Ra + hyperfine constants represent the best benchmark 
testing cases for Th 3+ , therefore, we carry out new cal- 
culations for Ba + and Ra + hyperfine constants A using 
exactly the same approach used for Th 3+ . The lowest- 
order values are listed in the "DF" column of Table Q] The 
relative correlation correction, evaluated as the difference 
of the ab initio final and the lowest-order values, is listed 
in column "CC" . We find excellent agreement with exper- 
imental values with exception of the nd 5 / 2 states, where 
the correlation correction is so large that the final and 
lowest-order values are of opposite sign. Our Fr result 
for the A(9s) state has been used in Rcf. [20] to extract 
the value of the 210 Fr nuclear magnetic dipole moment 
with 1% uncertainty. 

In this work, we have carried out the calculations of 
the hyperfine constants for the 5/, 6d, 7s, 7p, 7d, 8s, 
and 8p states using an SDpT all-order method. A Fermi 
distribution with the same parameters as the charged dis- 
tribution is used in all calculations to model a finite mag- 
netization distribution. The Breit interaction has been 
added on the same footing as the Coulomb interaction in 
the construction of the finite basis set. It was found to 
be small, less than 1%. 

Determination of 229 Th (I — 5/2) nuclear magnetic 
dipole moment is illustrated in Table [II] The values of 
fi (in (1b) are obtained as the ratio of the experimental 



TABLE II: Determination of the Th nuclear magnetic 
dipole moment. The values of u (in /is) are obtained as the 
ratios of the values in the Expt. and Aj\i theory columns. 
The value determined from the 6^3/2 hyperfine constant is 
taken as the final value. Theoretical values of A (in MHz) 
with /1 = 0.360 fiB are listed in column "A, Theory". Rel- 
ative differences of the experimental hyperfine constants and 
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A values and theoretical A/ n results for all four states. 
Ab initio theory values are listed for the 5/ 5 / 2 and 5/7/2 
states. The values for 6d 3 / 2 and 6g? 5 / 2 states are obtained 
by scaling ab initio results by -0.75% and -13.7%, respec- 
tively, using a comparison of theoretical bd Ba + values 
with experimental data. The scaling factors were calcu- 
lated as the difference of the value for the corresponding 
bd Ba + constant with experiment multiplied by the ratio 
of the correlation corrections for the 6dj state of Th 3+ 
and 5dj state of Ba + . Ba + is an excellent benchmark 
case for comparison with Th 3+ since the distribution of 
the correlation correction terms for all states is very sim- 
ilar for the two and the magnetic moment of 137 Ba is 
known with less than 0.01% uncertainty (22|. The value 
determined from the 6d 3 / 2 hyperfine constant is taken 
as the final result because of the Ba + benchmark com- 
parison that allowed for additional improvement of the 
value. We note that the value extracted from the 5/ 5 / 2 
state differs by only 0.38% from the final value. The 
5/ 5 / 2 ab initio result is expected to be the most accu- 
rate since it has the smallest correlation correction, only 
11%, while the correlation correction for the 6d 3 / 2 state 
is 24%. To illustrate the consistency of the extracted /i 
values with all four experimental results, we multiply the 
theory values for A//i by our final value of /i = 0.360/is 
and compare them with experimental data (6d3/ 2 will be 
in exact agrement since it was used to determine the final 
value of n). The relative differences of theory and exper- 
iment are listed in the last column "Diff." in %. We find 
that all four of the theory values are in agreement with 
experiment within experimental uncertainties listed in % 
in column "Unc." for convinience. Our final value of 
fi = 0.360(7) hb is significantly lower than the previous 
result fj, — 0.46(4) [is [2J] that was also derived from the 
hyperfine constant measurements. 

Determination of the 229 Th nuclear electric quadrupole 
moment is illustrated in Table fTTTl The table structure is 
similar to that of Table [TH but also includes prior deter- 
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TABLE III: Determination of the 229 Th nuclear electric quadrupole moment from theory values of Ref. and present results. 
The "Present work" B/Q column contains final ab initio theory results in MHz/eb. The average value of Q is taken as the 
final value. Theoretical values of B with Q = 3.11 eb are listed in columns "B" . The relative difference of the experimental 
hyperfine constants and theoretical values with Q — 3.11 eb are given in column "Diff." in %. 
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TABLE IV: Hyperfine constants A and B in 229 Th 3+ (1=5/2, 
fj, — 0.360 fj,B, Q = 3.11 eb) in MHz. Relative correla- 
tion corrections are listed in columns "CC" in %. '"'These 
three values are obtained by scaling ab initio 156.5, -11.4, and 
115.4 MHz results using Ba + data. 
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mination of Q carried out in .10] using theory data from 
[j| for comparison. Columns 3-6 show determination of 
Q using theory values from Q, and columns 7-10 show 
determination of Q using our present ab initio theory re- 
sults. We take the average of the results obtained from 
four states as the final Q value since the experimental 
and theoretical accuracy is similar for all four states and 
no relevant benchmarks exist for any of the four states. 
The value of the Ba + nuclear quadrupole moment [22[ 
relies on the old many-body perturbation theory calcu- 
lations which are of lower accuracy than the present cal- 
culation, so Ba + experimental B values can not be used 
as a benchmark. Remarkably, the average value of Q 
determined using our data and Ref. [|| is identical to 3 
significant figures, despite 2-6% differences in the values 
of B/Q. However, the differences of Q values determined 
from all four 5fj and 6dj states and the final value are 
lower in the present work, 0.7-1.3% (compare two "Diff". 
columns). Our final value, Q = 3.11(6) eb, is signifi- 
cantly lower than the previous values of 4.3(9) eb [13] 
also inferred from optical hyperfine measurements, but is 
in agreement with the result deduced from the Coulomb 
excitation of the nucleus, Q = 3.15(3) eb [35|. 

We used several different methods to establish the un- 
certainties in the values of \i and Q. 

(Ia) First, we establish the size of the total correlation 
corrections to the 5/ and 6d hyperfine constants. We 



list 229 Th 3+ hyperfine constants A and B in Table ITVl in 
MHz; our final values /i = 0.360 hb, Q = 3.11 eb are 
used. The relative size of the correlation corrections de- 
fined as the ratio of ab initio (SDpT-DF)/SDpT, where 
DF are the lowest-order results, is listed in columns "CC" 
in %. We find that the correlation corrections to the 6d 
and 5/ A and B values are relatively small, 11-29%, with 
the only exception of A(6d 5 / 2 ), where the correlation cor- 
rection is so large that the lowest-order DF value has a 
sign opposite that of the final value. The uncertainty of 
both theory and experiment is about 5% for the 6d 5 / 2 
state, so it was not used in determination of \i. 

(Ib) Next, we establish the clear correlation between 
the size of the correlation corrections and the accuracy 
of the values using the comparison of the theory and ex- 
perimental values for Ca+, Sr+, Hg + , Rb, Cs, Ba + , Fr, 
and Ra + , given in Table I of the supplementary material 
[h|. All cases similar to the 5/ and 6d states of Th 3+ 
where the correlation correction is below 30% agree with 
experiment to better than 2%. We note that Hg + is not 
a true "monovalent" case owing to low-lying core exci- 
tations that are not present in Th 3+ , so agreement is 
slightly worse (-2.3%). We find the average ratio of the 
last two columns of Table I of [15j (that give a correla- 
tion correction and difference with experiment) is 2.9%, 
and the maximum ratio is 6% [we excluded Hg + and 
the anomalous case of the 4d 3 / 2 state in Rb, where the 
correlation is 65% for the average]. If we average over 
only heavy systems listed in Table H](Cs, Ba + , Fr, Ra + ) 
the average ratio is still 3.1%. Therefore, we can make 
a general conclusion that the uncertainty of the calcula- 
tions is expected to be on the order of 3% of the total 
correlation correction and should not exceed 6%. For 
the 6d 3 / 2 and 5/j A constants, 3%/6% of the correla- 
tion gives (0. 3-0. 7%)/(0. 6-1.4%) uncertainty, depending 
on the state, with 5/5/2 being the most accurate. For all 
6dj and 5/j B constants, 3%/6% of the correlation gives 
(0.5-0.9%)/(l-1.8%) uncertainty, depending on the state, 
with 6g?3/2 being the most accurate. This analysis gives 
1.5% upper bound on the theory uncertainty for [i and 2% 
for Q. While the above analysis is done for A constants, 
it holds for a variety of states with completely different 
distributions of the correlation correction terms. There- 
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fore, it is expected to be applicable to B constants as 
well. The experimental uncertainty is 0.7% for A(6d 3 / 2 ) 
and .4(5/5/2), giving a 2% combined upper bound on the 
uncertainty of /x. The experimental uncertainty for all B 
constants is 0.26-0.47%, which can be considered negligi- 
ble in comparison with the theory uncertainty. 

(II) We further compare a benchmark case of A(6d 3 / 2 ) 
with A{nd 3 / 2 ) Ba + and Ra + measurements since our fi- 
nal determination of \x is based on this value. To ascer- 
tain that these cases are indeed very similar, we com- 
pare the entire breakdown of the 20 correlation correc- 
tion terms and normalization contributions to the hyper- 
fine constants. The comparative breakdown of correla- 
tion terms for various states of Ba + , Ra + , and Th 3+ is 
given in Table II of the supplementary material 15j . We 
find that the Ba+ and Ra + nd 3 / 2 cases are nearly iden- 
tical, while the Th 3+ 6d 3 / 2 case is very similar, but has 
a smaller overall correlation correction, 24%, instead of 
33-34%. This is expected as the correlation effects de- 
crease with increasing degree of ionization. As we noted 
above, we use a —1% difference of our Ba + value with 
experiment and a 24/34 correlation correction ratio to 
adjust our 6d 3 / 2 value by —0.75%. As a result, we ex- 
pect that our 6^3/2 theory value is accurate to better 
than 1%. However, the difference for the similarly ad- 
justed 6d 3 / 2 hyperfine constant of 213 Ra + with experi- 
ment is —1.9%. The Ra + measurement uncertainty is 
1.1% for this state, and the uncertainty of the 213 Ra + 
magnetic moment (from a direct measurement) is 0.3%. 
We note that there is a 1% inconsistency in the values of 
measured 7s hyperfine constants in 211 Ra + and 213 Ra + 
isotopes and our corrected (using Ba + data) values of the 
6P3/2 constant is in perfect agreement with the experi- 
ment. Therefore, it is likely that the discrepancy is due 
to uncertainty in the Ra + measurement. Nevertheless, 
we keep our original 2% estimated uncertainty of ji. 

(III) Finally, we use the consistency of /i and Q ob- 
tained from different states as an independent uncer- 
tainty estimate. Table II of the supplementary mate- 
rial [15| illustrates that the contributions of the various 
correlation terms are quite different for the 6d 3 /2, 5/ 5 / 2 , 
and 5/7/2 states. In fact, the correlation correction is 
negative for the 5/7/2 level while it is positive for the 
other two levels. Therefore, these calculations are suffi- 
ciently different that the spread of the fi values obtained 
from three levels (0.8-1.3%) provides another indepen- 
dent estimate of the accuracy. The 2% uncertainty on 
the 4(5/7/2) measurement does not presently allow us 
to significantly reduce the uncertainty. The comparative 
breakdown of the correlation correction to the B con- 
stants for the 6d and 5/ levels given in Table III of the 
supplementary material also shows differences between 
the correlation terms for the 6d and 5/ levels. The 0.7- 
1.4% difference of the Q values obtained for four levels 
with the average final value confirms a 2% uncertainty 
estimate for Q. 

The present accuracy of the value of /z can be improved 
(to 1% uncertainty or better) by remeasuring A for the 



6d 3 / 2 and 5/j levels to 0.3% or better. If [i derived from 
all three states were found to be consistent to 0.5-1% with 
improved experimental accuracy, this would demonstrate 
higher accuracy of the theory values and, subsequently, 
\x. The measurements of the 7s and 7p\/ 2 hyperfine con- 
stants A can be used to better understand finite mag- 
netization distribution effects, serve as an independent 
benchmark test case for Fr and Ra + measurements, and 
may provide improvement in the accuracy of /i. Our cal- 
culations give 3.1% and 1.1% finite magnetization distri- 
bution corrections for the 7s and 7p\/ 2 states which may 
be a substantial overestimation based on rough nuclear 
calculations [2(J. This seems to be confirmed by differ- 
ences in comparison of Cs/Fr and Ba + /Ra + values for 
these states as finite magnetization distribution effects 
are much smaller in Cs and Ba + . While the uncertainty 
of our ab initio calculations for A(7p 3 / 2 ) is about 3-4%, 
our scaled value should be accurate to about 1%. The 
magnetization distribution effects arc essentially zero for 
all states other than the 7s and 7p\/ 2 states. 

Improved measurements of the 6d 3 / 2 and 5fj hyper- 
fine intervals are critical for further improvements of the 
nuclear magnetic dipole moment value. Using an ultra- 
cold sample sympathetically cooled with, e.g., 232 Th 3+ , 
and sub-kHz linewidth lasers, the ~ 100 kHz wide optical 
transitions of interest could all be measured with inaccu- 
racies of < 10 kHz. In extracting the hyperfine intervals, 
this would reduce the current uncertainties in the A con- 
stants from ~ 10 -2 to ~ 10~ 5 , a level well below that of 
the companion calculations. Such measurements would 
also reduce the current uncertainties in the measured B 
constants by 3 orders of magnitude. 

The structure of the 229 Th 3+ nucleus may also be stud- 
ied beyond the first two electromagnetic moments. The 
5/5/2 electronic ground level is well suited for such stud- 
ies, as it has a large total angular momentum. Because J 
is greater than 1, the valence electron has non-zero cou- 
pling to the nuclear magnetic octupole moment n. This 
effect, parameterized by the magnetic octupole hyper- 
fine constant C, is expected to shift hyperfine intervals 
at the level of 10 1 — 10 2 Hz. Observing such an effect is 
straightforward in this ground-state hyperfine manifold, 
using microwave spectroscopy, as the individual states 
are effectively infinitely narrow and each hyperfine level 
contains an mp = clock state, allowing for the removal 
of first-order Zeeman shifts from the measurements. All 
hyperfine intervals should be measurable with inaccura- 
cies of ~ 0.01 Hz. This level of error is < 1% of an 
expected minimum value for C, indicating that percent- 
level inaccuracies in atomic structure calculations and 
nuclear magnetic dipole and electric quadrupole moment 
values would lead to extraction of n with an uncertainty 
of at most a few percent, on par with the most accurately 
determined nuclear octupole moment to date (36j. 

In conclusion, we determined the nuclear magnetic 
dipole and electric quadrupole moments of 229 Th with 2% 
uncertainty by combining high-precision theoretical and 
experimental values of the hyperfine constants and out- 
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line a pathway for further improvements in the accuracy 
of these properties. The present work also establishes a 
high-precision benchmark for accurate determination of 
electric quadrupole nuclear moments in other systems. 
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